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Superconductivity in Pb, H under extreme pressure and CaBeSi, in tlie framework of the density 
functional theory for superconductors, is discussed. A detailed analysis on how the electron-phonon 
and electron-electron interactions combine together to determine the superconducting gap and crit- 
ical temperature of these systems is presented. Pb, H under pressure and CaBeSi are multigap 
superconductors. We will address the question under which conditions does a system exhibits this 
phenomenon. The presented results contribute to the understanding of multiband and anisotropic 
superconductivity, which has received a lot of attention since the discovery of MgB2, and show how 
it is possible to describe the superconducting properties of real materials on a fully ab-imtio basis. 



I. INTRODUCTION 

The possibility of multiband superconductivity was proposed [l| a few years after the formulation of BCS theory 
and was later experimentally observed 0, [3, 13 • this regard, the details of the band (n) and k dependence of the 
superconducting (SC) gap were analyzed for many materials like Pb [5, 6, 7] and Sn Q. The recent discovery of two- 
gap superconductivity in MgB2 has renewed interest in this phenomenon. Theoretical and experimental investigations 
have emphasized the importance of multiple gaps in the enhancement of the critical temperature . Theoretical 
estimates have shown that including the information about the cr, tt band dependence of the electron-phonon (e-ph) 
interaction in MgB2 greatly enhances the calculated Tc and improves the agreement with experiment, as compared 
to an averaged, single-band, calculation]^ [l^, [ll|, [l^. Therefore of crucial importance is the question: When does 
the presence of multiple gaps occur and how does it affect Tel In this context, it is very interesting to investigate the 
different materials showing this peculiarity. 

Recently, a novel approach to superconductivity, based on density functional theory (SCDFT)[l^, has been 
able to describe, in a completely parameter-free fashion, the superconducting properties of several materials, ranging 
from the weak and intermediate to strong coupling regime p^. Il3l [T^ and from ambient to high pressure conditions 
[l^[l^. Unlike in the Eliashberg theory, in SCDFT the Coulomb interaction is treated on the same footing as the 
electron-phonon interaction. This is crucial for determining Tc accurately. 

In this work, we solve the SCDFT gap equation in a fully (n, k)-resolved formalism and calculate from first-principles 
the SC properties of Pb, H under pressure and CaBeSi . Within this approach all the features of the SC gap emerge 
naturally, without any a priori, material specific physical model or adjustable parameters. Our results confirm the 
experimental finding that Pb is a two band superconductor. The calculated gaps, their overall anisotropy, and the 
Tc result in good agreement with experiments reported in Refs. [13, [H- A Tc enhancement of ~ 8% can be directly 
related to the presence of multiple gaps f\^. 

The SCDFT, with its truly ab-initio character, allows not only to provide a better understanding of experimental 
facts, but also and more importantly, to predict the superconducting nature of systems which are not yet experimen- 
tally realized. One prime example among those is H under extreme pressure. The possibility of high-temperature 
superconductivity in metallic hydrogen was sugg ested in 1968 by N. W. Ashcroft [20]. Since then it has represented 
one of the most fascinating and intriguing topics involving fundamental issues ranging from the limits of phonon 
mediation in the superconducting phenomenon, to implications in astrophysics [21j. In the following we show that 
SCDFT predicts electron-phonon mediated superconductivity in H, with Tc values up to 242 K at 450 GPa. 

Another system for which SCDFT predicts multigap SC is CaBeSi [23|, a material isostructural isoelectronic to 
MgB2. Despite the many similarities with MgB2, according to our calculations CaBeSi has a very low critical 
temperature (Tc~ 0.4 K), consistent with experiment [1^. However, CaBeSi exhibits a complex gap structure, with 
three distinct gaps at the Fermi level. 

The paper is organized as follows: In Section [TT] we summarize the main features of SCDFT and describe our 
computational approach; in Sections IIIIl IIVI and |V] we present our results for Pb, H under pressure and CaBeSi 
respectively; finally, in Section [VII we summarize our conclusions. 
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II. THE DENSITY FUNCTIONAL THEORY FOR SUPERCONDUCTORS 

Density functional theory (23 | has enjoyed increasing popularity as a reliable and relatively inexpensive tool to 
describe real materials. In this section we will briefly outline the DFT approach to superconductivity, and refer 
to the original papers for more details. In order to give an introduction to SCDFT, it is instructive to recall how 
magnetism is treated within DFT. The Hohenberg-Kohn (HK) theorem [2^ states that all observables, in particular 
also the magnetization, are functionals of the electronic density alone. This, however, assumes the knowledge of the 
magnetization as a functional of the density. Finding an approximation for this functional is extremely hard and, in 
practice, one chooses a different approach. The task can be vastly simplified by treating the magnetization density 
m(r), i.e., the order parameter of the magnetic state, as an additional fundamental density in the density functional 
framework [l^. An auxiliary field - here a magnetic field Bext(r) - is introduced, which couples to m(r) and breaks 
the corresponding (rotational) symmetry of the Hamiltonian. This field drives the system into the ordered state. If 
the system is actually magnetic, the order parameter will survive when the auxiliary perturbation is quenched. In 
this way, the ground-state magnetization density is determined by minimizing the total energy functional (free energy 
functional for finite temperature calculations) with respect to both the normal density and the magnetization density. 
Within this approach much simpler approximations to the xc functional (now a functional of two densities) can lead 
to satisfactory results. 

The same idea is also at the heart of density functional theory for superconductors, as formulated by Oliveira, Gross 
and Kohn ^2f\. Here the order parameter is the so-called anomalous density, 

x(r,r') = (4'T(r)*,(r')>, (1) 

and the corresponding potential is the non-local pairing potential A(r, r'). It can be interpreted as an external pairing 
field, induced by an adjacent superconductor via the proximity effect. Again, this external field only acts to break 
the symmetry (here the gauge symmetry) of the system, and is quenched at the end of the calculation. As in the case 
of magnetism, if the system is actually a superconductor the order parameter will be sustained by the self-consistent 
effective pairing field. The approach outlined so far captures, in principle, all the electronic degrees of freedom. To 
describe conventional phonon-mediated superconductors, also the electron-phonon interaction has to be taken into 
account. 

In order to treat both weak and strong electron-phonon coupling, the electronic and the nuclear degrees of freedom 
have to be treated on equal footing. This can be achieved by a multi-component DFT, based on both the electronic 
density and the diagonal of the nuclear N-body density matrix [H, [2^ [s^ : 



r(R) = ($t(Ri) . . . $t(R^)$(R^) . . . $(R^))^ (2) 

where $(R) is a nuclear field operator. The quantity r(R) is then included as an additional "density" in the formalism, 
besides the electronic density and the SC order parameter. 

In order to formulate a Hohenberg-Kohn theorem for this system, we introduce a set of three potentials, which 
couple to the three densities described above. Since the electron-nuclear interaction, which in conventional DFT 
constitutes the external potential, is treated explicitly in this formalism, it is not part of the external potential. The 
nuclear Coulomb interaction [/"" already has the form of an external many-body potential, coupling to r(R), and 
for the sake of the Hohenberg-Kohn theorem, this potential will be allowed to take the form of an arbitrary N-body 
potential. All three external potentials are merely mathematical devices, required to formulate a Hohenberg-Kohn 
theorem. At the end of the derivation, the external electronic and pairing potentials will be set to zero while the 
external nuclear many-body potential is set to the nuclear Coulomb interaction. 

As usual, the Hohenberg-Kohn theorem guarantees a one-to-one mapping between the set of the densities 
{n(r), x(r, r'), r(R)} in thermal equilibrium and the set of their conjugate potentials {w°xt(r) — M, Aoxt(r, r'), w^xt (R)}- 
Therefore all the observables are functionals of the set of densities. Finally, it assures that the grand canonical po- 
tential, 

n[n,x,r]=F[n,x,T] + J d^r nirM^,{r) - f,] - J d\J d^r' [x(r, r')A:,t(r, r') + h.c] + J d^RTimv^.M), (3) 

is minimized by the equilibrium densities. We use the notation A[f] to denote that A is a functional of /. The 
functional _F[7i, x,r] is universal in the sense that it does not depend on the external potentials. It is defined by 

F[n, X, r] = T-[n, x, T] + T"[7i, x, T] + C/™[n, x, F] + C/-[n, x, F] - ^S[n, x, F] , (4) 
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where T*' represents the electronic kinetic energy, C/*^" the electron-electron interaction, T" the nuclear kinetic energy, 
[/"" the Coulomb repulsion between nuclei, the electron-nuclei interaction and S is the entropy of the system, 

S[n,x,T] -Tr{/5o[n,x,r]ln(/5o[n,x,r])}. (5) 

In standard DFT one normally defines a Kohn-Sham system, i.e., a non-interacting system chosen such that it has 
the same ground-state density as the interacting one. The variational procedure for this system gives Schrodinger-like 
(Kohn-Sham) equations for non-interacting electrons subject to an effective (Kohn-Sham) potential. These equations 
are nowadays routinely solved by solid state theorists. In our formalism, the Kohn-Sham system consists of non- 
interacting (superconducting) electrons, and interacting nuclei. We will not describe here the details of the method, 
and will only outline its basic features: The Kohn-Sham potentials, which are derived in analogy to normal DFT, 
include the external fields, Hartree, and exchange-correlation terms. The latter account for all many-body effects of the 
electron-electron and electron-nuclear interactions. Obtaining their explicit form has represented a major theoretical 
effort [30I Isil . [33 |. Once this problem has been solved, the problem of minimizing the Kohn-Sham grand canonical 
potential can be transformed into a set of three differential equations that have to be solved self-consistently: One 
equation for the nuclei, which resembles the familiar nuclear Born-Oppenheimer equation, and two coupled equations 
which describe the electronic degrees of freedom and have the algebraic structure of the Bogoliubov-de Gennes [33l | 
equations. 

The resulting Kohn-Sham Bogoliubov-de Gennes (KS-BdG) equations read (we use atomic Hartree units) 



" 2 
2 



+ <(i-)-M 

+ <(r) - At 



d^r' As(r, r')w„k(r') = ^„k w„k(r) , 
dV A:(r,r')w„k(r') = ^„kW„k(r) , 



(6a) 
(6b) 



where w„k(r) and w„k(r) are the particle and hole amplitudes. This equation is very similar to the Kohn-Sham 
equations in the OGK formalism [23|. However, in the present formulation the lattice potential is not considered 
an external potential but enters via the electron-ion Hartree term. Furthermore, our exchange-correlation potentials 
depend on the nuclear density matrix, and therefore on the phonons. Although Eq.Q and the corresponding equation 
for the nuclei have the structure of static mean-field equations, they contain, in principle, all correlation and retardation 
effects through the exchange-correlation potentials. 

These KS-BdG equations can be simplified by the so-called decoupling approximation [l^ [s^l , which corresponds 
to approximating the particle and hole amplitudes by: 

"nk(l") ~ Wnk<y5„k(l") ; W„k(r) ~ Vnl^^ril^{r) , (7) 

where the wave functions ifuki^c) are the solutions of the normal Schrodinger equation. In this way the eigenvalues in 
Eqs. ([6]) become Enk = i-^nk, where 



nk 



(8) 



and ^„k = ^nk—fJ-- This form of the cigcnenergies allows us to interpret the pair potential A„k as the gap function of the 
superconductor. Furthermore, the coefficients u„k and u„k are given by simple expressions within this approximation 



Unk = --^sgn(i5„k)e"^ 



'1 



Vnk 



711 



Enk 



Finally, the matrix elements A„k are defined as 

A„k = Jd\J dV <kWAs(r,r')^„k(r'), 
and (pnk is the phase e"^"'' = A„k/|A„k|- The normal and the anomalous densities can then be easily obtained from: 



(9a) 
(9b) 

(10) 
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(11a) 
(lib) 
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FIG. 1: Superconducting gap of Pb as a function of the energy distance from the Fermi energy (T = K). SmaU black points: 
A resulting from the inclusion of the k-resolved e-ph and Coulomb matrix elements. Small orange points: A calculated with 
a corresponding averaged e-ph interaction (Eliashberg function a^FlQ,)). Big triangles and squares: experimental values (from 
Refs. 113,1181). 



Within the dccouphng approximation, we finally arrive at an equation for the k-resolved superconducting gap A„k) 
which has the foUowing form[ll [ll [23, [Hi: 

1 tanh (^§En'k'^ 

n'K' 

Eq. (fT2)) is the central equation of the DFT for superconductors. The kernel K. consists of two contributions 
K, — /C°^P'' + JC^~'^, representing the effects of the e-ph and of the e-e interactions, respectively. The diagonal term 
Z plays a similar role as the renormalization term in the Eliashberg equations. Explicit expressions of JC~v^ and Z, 
which are the results of the approximate functionals, are given in Eqs. 9 and 11 of Ref. respectively. These two 
terms involve the e-ph coupling matrix, while 1C'^~'^ contains the matrix elements of the screened Coulomb interaction. 
Eq. (|12p has the same structure as the BCS gap equation, with the kernel JC replacing the model interaction of BCS 
theory. This similarity allows us to interpret the kernel as an effective interaction responsible for the binding of the 
Cooper pairs. Moreover, we emphasize that Eq. p2|) is not a mean-field equation (as in BCS theory), since it contains 
correlation effects via the SC exchange-correlation functional entering K, and Z. Furthermore, it has the form of a 
static equation - i.e., it does not depend explicitly on the frequency - and therefore has a simpler (and computationally 
more manageable) structure than the Eliashberg equations. However, this certainly does not imply that retardation 
effects are absent from the theory. Once again, retardation effects enter through the xc functional, as explained in 
Refs. [13, 14] . The SCDFT allows to treat the e-ph and the screened e-e interactions on the same footing. These 
terms, however, can be treated at different levels of approximation. 

We calculated the screened Coulomb matrix elements (ME) with respect to the Bloch functions, for the whole 
energy range of relevant valence and conduction states. The different nature of the electronic bands in each material 
{e.g. some of them can be highly localized while others more delocalized) , strongly calls for the use of a non-diagonal 
screening, including local field effects. 

As mentioned above, the normal state calculations, necessary for the study of the superconducting state, are 
performed within DFT in the LDA or GGA approximations. Computationally, the electronic and phononic properties 
are obtained using the pseudopotential method as implemented in the QUANTUM-ESPRESSO package[36i|; the 
screened Coulomb matrix elements are obtained with the SELF[33 code. 



III. MULTIBAND SUPERCONDUCTIVITY IN Pb 



The peculiar aspects of superconductivity in Pb have been discussed in the past d, H, [3, [13, [H, [H, [H, [i3| ■ In 
this section we present the results for the SC properties of Pb, calculated within the SCDFT formalism and highlight 
the importance of the (n, k) gap anisotropy in this system. We first concentrate on the SC gap function, i.e. the 
solution of Eq. p^ . There are two ways to visualize this function: The first is to look at the gap as a function of 
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the normal-state energy eigenvalue. The set of black points in Fig. [T] shows the gap at T = K calculated on a set of 
random k points, as a function of the energy distance from the Fermi level Ep. An important feature of this plot is 
that for each energy the gap is not a single-valued function. This means that, in general, the SC gap is not isotropic 
in reciprocal space, i.e. for e =Ep its value depends on the Fermi vector k^^. In particular, for each e ~ Ep we 
observe two distinct "sets" of gaps, in strict analogy with the case of MgB2 (see, for comparison, Fig.l in Ref. 
Moreover, each "set" has an associated, finite, vertical energy spread. 

Most importantly the gap function has a negative tail extending to high energy. This is a general feature of the gap 
function, due to the presence of the Coulomb repulsion and to its dominance for large energies. Due to this feature. 
Coulomb renormalization effects are properly included in the SCDFT formalism. We emphasize that these effects are 
crucial to obtain a predictive superconducting solution. 

A second possibility to visualize the gap is to define equi-energy surfaces and plot the gap value on that surface. 
The equi-energy surface can be chosen to be the Fermi surface (FS). In Fig. [2^ we plot the FS of Pb, consisting of two 
separate sheets, coming from the two (essentially p) bands crossing Ep. The sheets are topologically quite different, 
the first one being rather spherical and the second having a more complex tubular-like structure. The colorcrin Fig. 
[2^ represent the values of Ankp- It is clear that the two distinct sets of A in Fig. [l]come from the two sheets of the 
Fermi surface. Fig. [1] shows that their energy separation is in good agreement with the experiments [TtI. [isj. 

Fig. [5^ also shows that A„fc is anisotropic inside each single sheet, resulting in the vertical energy spread mentioned 
before (intraband gap anisotropy). Continuing our analysis, in Fig. we plot the e-ph coupling Xnkp on the FS, 
where 

I rm' 12 

A„fc = 2 V X Sisn'w - Ep) (13) 

777 i^k'-k.u 
n W M 

is the average of all the possible e-ph scattering processes connecting two points at FS, but always involving the 
electronic initial state (n, k). The total e-ph coupling constant A can be expressed as A = n(Ef) Snk -^nkX ^(cnk— ^'f)- 
Comparing the plots in Fig. [2^ and Fig. ^Bp we notice a striking similarity, that indicates a strong correlation between 
the ri-dependence and k anisotropy of the e-ph coupling and of the gap. This is a general feature that we also observe 
for other materials. Moreover, our calculations show that the gap anisotropy is due to the e-ph coupling only: the 
open orange symbols in Fig. [T] shows the gap function calculated with an average e-ph interaction, i.e. including the 
phononic kernels /C°~p'^ and Z that contain the Eliashberg function a'^F{VL) (see Eqs. 23 and 24 of Ref. [Til), but 
still with the Coulomb ME clear that the e-ph average washes out all the band and k dependence of 

the gap (i.e. A„fc = A(e„fc)). 

Eq. allows to calculate A„fc as a function of temperature T, defining Tc from the condition that /^nk{Tc) — 0. 
We obtain Tc= 5.25 K and Tc= 4.84 K for the anisotropic and average calculation respectively. We see that, although 
certainly not in a proportion comparable with MgB2, the presence of an anisotropic, multiband gap produces an 8% 
enhancement of the value of Tc . In MgB2 this effect is much more pronounced 0, El, than in the present case. 
This is partially due to the much lower difference in e-ph coupling between the two bands (in MgB2 the coupling in 
the a bands is roughly 3 times higher than in the tt bands). 
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FIG. 3: Fermi surface at 414 GPa. Different colors represents different values of the superconducting gap on the FS. Red 
(18.0 < A < 21 meV), yellow (13.6 < A < 18.0 meV) and blue (10.0 < A < 13.6 meV). 

An interesting question is under which conditions does a system show multiband superconductivity, i.e. with a 
clear separation between the gaps. Looking again at the Xnkp we see that in Pb the two sets corresponding to the 
two bands crossing Ep are not continuous, i.e. the values of Aifc^ and X2kF do not superimpose. In order to check 
whether the presence of two separate gaps is related directly to this property of the Xnkp , we performed several model 
calculations, rescaling the e-ph ME, in order to have a continuous set of Xnkp and keeping constant the value of the 
total, average A. Most interestingly we found that the two SC gaps Ankp get continuous exactly when the two set of 
Xnkp do, independent of the details of how the coupling were modified to bring the Xnkp to merge. This leads us to 
conclude that Pb is a two gap system due to the fact that Xnkp are disjoint sets (relative to the band n). In general 
there are two necessary conditions to fulfill in order to obtain multigap superconductivity. First the system needs to 
have disconnected FS sheets and second the electron-phonon coupling must be different on each of these sheets. 

IV. SUPERCONDUCTIVITY IN H UNDER PRESSURE 

Based on the simple BCS theory of superconductivity, we can understand why molecular metallic hydrogen could 
be a good superconductor |4l|: This system has very high phonon frequencies due to the light H mass, and it has 
a possibly strong electron-phonon interaction related to the lack of core electrons and to the quite strong covalent 
bonding within the H2 molecules. Many studies, aiming at investigating further this possibility [i^, E^, IS] , collected 
strong evidence pointing to high-Tc superconductivity. At present, however, the full scenario is still far from being 
clear and well established. 

The low temperature and high pressure (> 400 GPa) phase of molecular hydrogen is predicted to be a base- 
centered orthorhombic metallic molecular solid (known as Cmca phase ^fl.S]) with two molecules per unit cell located 
on different layers. The electronic band structure arises from the bonding and anti-bonding combination of the H2 
molecular orbitals. At high pressure, the band overlap between the valence and conduction bands produces a rich 
and complex Fermi surface (Fig. [3]) with disconnected sheets of different orbital nature[46] and provides the strong 
e-ph coupling necessary to superconductivity (see below) . 

Three types of bands form the Fermi surface of H at high pressure: The two tubular structures intersecting the ky 
axis and the two prism-like structures along the kx axis have hole-type character; the disk centered at the F-point and 
the structures near the C-point have electron-type character; the remaining sheets are related to degenerate orbitals 
along Z - T. The different orbital character of the FS branches suggests the occurrence of different couplings between 
the various bands, leading to an anisotropic SC gap and to multiband superconductivity. 

The molecular nature of the charge density distribution gives rise to an intriguing phonon dispersion (Fig[?]). The 
phonon modes can be grouped into three main branches: phononic, libronic and vibronic, corresponding respectively 
to the relative translations, rotations, and internal vibrations of the H2 molecules. The three FS regions possess 
distinct e-ph coupling: The most coupled one is the disk at F with 1.8 < Ak < 2.00; then, we find the two prism-like 
structures with 1.00 < Ak < 1.80, and the other sheets with Ak < 1-00. The tubular structures have the smallest 
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FIG. 4: H phonon dispersion at 414 GPa. 

Ak- The presence of multiple Fermi surfaces provides a " q-distributed" coupling, with many modes contributing to 
the pairing. This allows to increase A still avoiding a lattice instability, which could result from very high coupling 
at some specific q's. In fact, the Eliashberg function (FigIS]) shows three main frequency regions all of them strongly 
coupled to the electrons and associated with phononic, libronic and vibronic modes, in an increasing frequency order. 
However, a net distinction between phononic and libronic modes (in the low frequency region: < lu < 1600 cm~^) is 
prevented due to the branch mixing caused by the large renormalization of the libronic modes. As the pressure raises 
from 414 GPa, the coupling increases and a new feature appears at w 1700 cm^^ (Fig- El plot at 462 GPa). The 
additional peak in the a^F{uj) comes from an extra band which crosses the Fermi level close to the R-point. This 
creates a new electron-like Fermi sheet, strongly coupled with in-phase libronic modes, and opens a further coupling 
channel at higher pressure giving rise to a jump of A at w 435 GPa (inset of Fig. [S]). 

Although rather large, these values of A are of the same order as, e.g., those obtained in fcc-Li under high pressure 
(40 GPa). There, Tc reached ~ 20 K, a very large value for a free electron-like metal. In the case of hydrogen, we 
expect that the large phonon frequencies provided by the H2 vibronic modes lead to much higher T^s. 

The Tc was calculated solving the SCDFT anisotropic gap equation Eq. p^ . As previously shown in the case of 
low-density metals [l^ and discussed by Richardson and Ashcroftfl^] a complete inclusion of both electron-phonon 
and electron-electron interactions is crucial to reproduce the experimental T^s. We find very high T^s, considerably 
increasing with pressure up to 242 K at 450 GPa (inset of Figl6]). Comparing Figs. [5] and [6] insets, it is clear that Tc 
follows closely the A behavior with pressure. 

We predict a multigap (three gap) superconductivity, which is progressively lost at high pressure, where interband 
scattering raises the nondiagonal matrix elements, leading to a merging of the different gaps. The occurrence of 
multigap superconductivity is clearly illustrated in Fig. [SI which shows the A„fc at the Fermi level as a function of 
temperature [for all the bands (n) and k-points] and at the pressure with largest band anisotropy, P — 414 GPa (Tc 
= 84 K). The three gaps are evident: The largest one at « 19.3 meV, and two very anisotropic, overlapping gaps 
at 15.4 and 13.6 meV. The three gaps are associated with the three different FS sheets (cf. Fig. [3]): The strongly 
coupled disk around the F-point, the prism-like sheets, and others associated with the lowest gap. 
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FIG. 5: a^F(f2) at 414 (solid line) and 462 (dashed line) GPa. The inset shows the electron-phonon coupling (A) as a function 
of the pressure. 
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FIG. 6: A,jk at Ef at 414 GPa. The shaded regions represent the k-anisotropy over different bands (note that the blue 
and yellow gaps overlap from 13.00 to 15.75 meV). The inset shows the superconducting critical temperature as a function of 
pressure. 




FIG. 7: CaBeSi Fermi surface. The colors identify the three gaps at Ep: the largest a (yellow and blue cylindrical-like sheets); 
the intermediate n bonding (green spheres); the lowest n antibonding (violet). 

In H the main source of band anisotropy is given by the e-ph interaction. Furthermore a simple multigap BCS 
model (not including anisotropy in k) although not adequate to treat strong coupling system, is able to reproduce 
qualitatively the full SCDFT results [40]. This suggests that the dominant effect of the anisotropy on Tc is related 
to the band contribution and that the k space anisotropy is less relevant. We emphasize, however, that although 
the anisotropy of the interaction generates three different gaps, it does not play a fundamental role in enhancing 
Tc, which shows an increase of only lOK compared to the value obtained within an isotropic approximation. This 
is due to the fact that in H the inter-band (low-q) coupling dominates over the intra-band coupling (high-q) [48| . 
This is different from MgB2, where the intra-band coupling is much stronger and the anisotropy doubles the critical 
temperature [Tol. [T^. 

V. SUPERCONDUCTIVITY IN CaBeSi 

In this section we report the SC properties of CaBej.Si2_a; at a; = 1 (CaBeSi in the following) in the AIB2 phase, 
with Be and Si atoms alternating within the honeycomb layers with an AA stacking. A general similarity is found 
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FIG. 8; CaBeSi superconducting gap as a function of the energy distance from the Fermi leveL 

between CaBeSi and MgB2 : in both materials the ct, 7r-bonding (tt;,) and 7r-antibonding (77^) bands cross the Fermi 
level (Ep). In CaBeSi the nb bands are almost fully occupied, leaving only small hole pockets at K which give rise to 
the little tt spheres of the Fermi surface (Fig. [7]). These replace the tt tubular structure present in MgB2. 

Despite the general similarities, CaBeSi and MgB2 have rather different chemical and phononic properties, with 
important consequences for the e-ph interaction [22|- In fact, CaBeSi phonon frequencies are lower in comparison with 
MgB2, mainly due to the larger mass of Ca and Si versus Mg and B. In the CaBeSi phonon dispersion, the E2g mode 
is fairly flat along the in-p lane BZ symmetry lines and it shows only a very weak renormalization (four times smaller 
than in MgB2 [491 [sol. fsifl due to the very small E2g electron-phonon matrix elements. In turn, this is related to the 
delocalized and ionic nature of the a bonds in CaBeSi . In fact, the connection between strongly covalent bonds and 
strong e-ph coupling seems to be a general feature [3, [Tg, 19, 52]. The small e-ph coupling makes CaBeSi a weak 
coupling superconductor with e-ph A = 0.38 and Tc— 0.4 K, in spite of a nesting at the Fermi surface twice as big as 
in MgB2 . 

The solution of the self-consistent gap equation reveals an unexpected complex structure (Fig. [5]) with clearly 
separated three gaps at Ep (unlike in H where two of the three gaps superpose). The calculated critical temperature 
is very low (Tc= 0.4), lower than the upper limit (4.2 K) set by the exp erimental results[2^. Unlike in MgB2, 
where superconductivity is interpreted within a two-band modeljl, IH, |54| . [sst . in CaBeSi there is a further nt-TTa 
gap splitting. As in MgB2, the largest gap is related to the cr FS sheets (cylindrical- like structures in Fig. [7]), the 
intermediate one to irt sheets (small hole spheres) and the lowest to TTa sheets. The additional nb-TTa gap splitting 
is a peculiar feature of CaBeSi not present in MgB2, where the two a and the two tt gaps merge together. In 
order to understand the origin of this splitting we perform some additional computational experiments, solving the 
gap equation (i) completely neglecting the Coulomb interaction, (ii) including only an averaged Coulomb term and 
[iii) with isotropically averaged Coulomb and phononic interactions, corresponding to the dirty limit. In both (i) 
and {ii) cases, the three gap structure is destroyed, bringing back to a two-band, MgB2-like gap structure. In case 
[Hi), instead, superconductivity is completely lost. As a result, we predict superconductivity in CaBeSi only if the 
anisotropic structure of the interactions is included. While, as in MgB2, the ct-tt gap splitting is related to the different 
e-ph coupling in these bands, the further Tib-T^a splitting is a pure effect of the complex structure of the anisotropic 
Coulomb repulsion, acting in a different way on the -Kb and tTq states. 

VI. SUMMARY 

Three multi-gap superconductors were investigated in this work: Pb, H under pressure, and CaBeSi. Without 
assuming any ad-hoc model or adjustable parameters, we find that Pb is a two-gap superconductor, in agreement 
with available experiments and with previous theoretical results based on other methods. We show that the n and 
k gap anisotropy and the separation of the two gaps correlate strongly with the anisotropy of the electron-phonon 
interaction on the different sheets of the Fermi surface. This multi-gap character produces an enhancement of Tc 
relative to the isotropic case. 

For H, we find a strong increase of the critical temperature with increasing pressure. As the corresponding pressure 
regime is not yet experimentally accesssible our calculation of the astonishingly high critical temperatures represents 
a genuine prediction. The analysis of the electronic and phononic properties points to the origins of the strong en- 
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hancement of superconductivity with rising presssure: In particular: i) wc confirm a strong clectron-phonon coupling, 
increasing with pressure; ii) the presence of multiple Fermi surfaces with different character provides both strong 
intra-band (low-q) and inter-band (high-q) electron-phonon scattering. In this respect, in) a major role is played 
by the molecular rotational (libronic) and vibrational (vibronic) modes, which are strongly coupled with the inter- 
molecular charge. Finally, iv) we predict three superconducting gaps at the Fermi level, associated with Fermi surface 
branches having different character. 

We have also investigated the superconducting properties of CaBeSi, a material isostructural and isoelectronic to 
MgB2. While the band structures present strong similarities, the phonon structure and the e-ph interaction of the two 
systems differ substantially. In particular, the less localized a charge of CaBeSi leads to a dramatic reduction in the 
electron-phonon coupling, with a consequent reduction of the phonon renormalization seen in MgB2. Interestingly, 
CaBeSi exhibits three superconducting gaps at the Fermi level. The iTb-TTa splitting is a pure effect of the complex 
structure of the anisotropic Coulomb repulsion, acting differently on the Tr^ and tTq states. 

The role of multigap superconductivity in enhancing the critical temperature of superconducting materials has been 
known for a long time. However, the crucial challenge is to predict under which circumstances and for which materials 
the phenomenon of multi-gap superconductivity actually occurs. With the SCDFT approach which features a fully 
anisotropic description of the electron-phonon coupling and of the screened Coulomb repulsion on an ab-initio basis, 
this goal has been achieved. 
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